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g i ABSTRACT 

We present a multi-wavelength observational study towards the high-mass young 



stellar object G8. 68-0. 37. A single massive gas-and-dust core is observed in the 
(sub) millimeter continuum and molecular line emissions. We fitted the spectral 



energy distribution (SED) from the dust continuum emission. The best-fit SED 
suggests the presence of two components with temperature of T^ - 20 K and 120 

KS ■ K, respectively. The core has a total mass of up to 1.5 x 10 3 M Q and bolometric lu- 

b ■ 

minosity of 2.3 x 10 4 L Q . Both the mass and luminosity are dominated by the cold 

component ( T& — 20 K) . The molecular lines of C 18 O, C 34 S, DCN, and thermally ex- 
cited CH3OH are detected in this core. Prominent infall signatures are observed in 
the 12 CO (1-0) and (2-1). We estimated an infall velocity of 0.45 km s _1 and mass 
infall rate of 7 x 10~ 4 M Q year -1 . From the molecular lines, we have found a high 
DCN abundance and relative abundance ratio to HCN. The overabundant DCN 
may originate from a significant deuteration in the previous cold pre-protostellar 
phase. And the DCN should now be rapidly sublimated from the grain mantles to 
maintain the overabundance in the gas phase. 

Key words: stars: pre-main sequence — ISM: molecules — ISM: kinematics and 
dynamics — ISM: individual (G8. 68-0.37) — stars: formation 
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1 INTRODUCTION 

Gravitational infall, or core collapse can take place in high mass young stellar object s (YSOs) at 



early stages and continue all the way to the stage of Ultra Compact (UC) Hn regions 



Sollins et al. 



e 



2005) . As shown by theoretical works Tiji na & Adamsll996t 



Keto 2003; 



Yorke & Sonnhalter2002; 



Gong & Ostrikell2009L etc.), the infall motion is critical for initiating the high mass star forma 
tion and maintaining the accretion flow to feed the stellar mass during the subsequent evolu- 
tionary stages. However, further observations are still needed to better constrain the physical 
properties of the infall, including its spatial distribution, mass infall rate, chemical effect, and 
to understand its relation with other dynamical processes, including outflow, di sk accretion, 



and core fragmentation. In the recent decade 



2003 



Fuller et al 



2005; 



Wyrowski et al 



2006; 



extensive spectroscopic surveys 



Purcell et al. 



2006; 



Wu et al 



(e.g. 



Wu & Evans 



2007 ) have been per- 



formed towards the potential high mass YSOs throughout the Milky Way. As a result, many in- 
fall candidates have been identified based on their spectral signatures. These sources can serve 
as good candidates to study the massive star birth and gas dynamics in the molecular cores. 
In the mean time, strong outflows are also widely detected towards those massive cores (e.g. 



Beuther et al. 



2002; 



. The infall and outflow motions should be 



closely related and interacting with each other throughout the star formation history. 



G8. 68-0.37 (G8.68 h ere after) is a young high-mass star forming region at a distance of 4.5 



kpc IMueller et al 



covered by 



200 2]) . In this r egion, compact multiple gas-and-dust clumps has been dis- 



Longmore et al 



GHz methanol masers (Wa lsh et al 



2011, L 



1 here after). The dusty core is associated with strong 6.7 



1998), but has no radio continuum emission, indicating 



that high mass stars are already formed, but have not yet ionized its surrounding gas. Lll also 
detected a bi-polar outflow in CO (2-1). The outflow may be responsible for the shock inter- 
action traced by the extended 4.5 urn emission (Figure 2 ther ein). In the mean time, the ob- 



servation in HCO + (1 - 0) suggests a plausible infall motion JPurcell et al. 



2006) which should 



be examined quantitatively. To improve the understanding in physical and chemical proper- 
ties of this source, we performed a multi -wavelength study using both the single dish antennae 
and the interferometers. The next section introduces the observations and data reduction, Sec- 
tion 3 presents the general observational results. Section 4 describes the dust continuum and 
molecular line emissions, wherein the infall signature is specifically described in Section 4.3. A 
summary is given in Section 5. 
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Infall motion in G8.68-0.37 3 
2 OBSERVATIONS AND DATA REDUCTION 

2. 1 The single dishes 

In Figure 1 we show the central positions and the beam sizes of all the observations. We have 
observed G8.68 using three different single-dish telescopes. In 2005, we observed HCN (3-2) 
and H 13 CO + (3 - 2) from the James Clerk Maxwell Telescopej ( JCMT). The poin ting center was 
adopted to be the coordinate of the strongest methanol maser IWalshet all 19981. with a position 
accuracy of 1.8"), which is close to the continuum emission peak of LI 1 (cross in Figure 1). The 
mapping step is 10" (corresponding to 1/2 beam size), as shown in Figure 2b. 

In November 2009, we observed / = 1 - line of 12 CO, 13 CO, and C 18 using the 13.7 m 
telescope at the Purple Mountain Observatory 2 ] (PMO). The PMO observation contains a grid- 
mapping with a coverage of several arc minutes over the region of G8.68. In this paper, we only 
use the spectra at one point which is closest to the continuum peak, as shown in Figure 1. The 
PMO beam is much larger than the CSO and JCMT, and is significantly deviated from the con- 
tinuum peak. Nevertheless, the beam has well covered the emission regions of the continuum 
and molecular lines. The data can thus be used to trace the gas motion on a larger scale near 
the core. 

In May 2011, the 7 = 2-1 lines of the three CO isotopologues were observed from the 10 
m telescope at the Caltech Submiilimeter Observatory 3 ] (CSO). The 12 CO (2 - 1) is observed at 
five symmetric points around the continuum peak with an offset of +23" in the R.A. and Dec. 
directions. Their positions are shown in Figure 2a. All the single-dish spectra are discussed in 
detailed in Section 3.2.2. 

In Table 1, we present the basic observational parameters and weather conditions for the 
three instruments. In Table 2, we shows the more specific observational parameters for the 
molecular lines. All the observations were performed in good weather conditions, with pointing 
accuracies better than 5". GILDAS software packaged is used for the data reduction and image 
plot. 

To measure the flux densities at different wavelengths, we also retrieved the Spitzer archival 

fl 
images at four IRAC bands from the GLIMPSE surveys, and the 24 and 70 urn images from 

1 JCMT is operated by the JAC, Hawaii, on behalf of the UK PPARC, the Netherlands OSR, and the Canadian NRC, see 

http://www.jach.hawaii.edu/JCMT/ 

2 http://www.dlh.pmo.cas.cn/ 
http://www.submm.caltech.edu/cso/ 

4 http://iram.fr/IRAMFR/GILDAS/ 

5 Availableathttp://irsa.ipac.caltech.edu/ seealso 



Benjamin et al 



2003) 
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and JCMT/ SCUBA images at 450 and 850 urn bands, which are available at the 



Canadian Astronomy Data Center (CADC) repository of the SCUBA Legacy Fundamental Ob- 
ject Catalogue^. 



2.2 The Submillimeter Array 

The Submillimeter Arraya (SMA) observations towards G8.68 are taken from the released SMA 
data archive. The observations are made in three epochs in the year 2007, 2008, and 2009, re- 
spectively. The observational parameters, including the calibration sources for each epoch are 
presented in Table 3. In all three epochs, the compact array was used, and the phase tracking 
center is RA.(J2000)=18 h 06 m 23.23 s , Dec.(J2000)=-21°37'14.19". The three observations have 
similar beam sizes for the synthesized and primary beams. In Figure 1 , we only show the beams 
of the 2008 observation in order to have a clear appearance. The observed gas-and-dust struc- 
tures (Figure 4) turn out to be smaller than the SMA primary beam. Thus the beam-edge weak- 
ening is not significant. The calibration and imaging were performed in Miriad 1 . The absolute 
flux level has an uncertainty of ~ 15%. The continuum emission was subtracted from the line- 
free channels in each sideband. The gain solution is self- calibrated for the continuum image 
and then exported to the spectral line data. 

We note that among the SMA data, the 2008 observation has the longest on-source integra- 
tion time hence the lowest noise level. In addition, in 2008 all eight antennae of the SMA were at 
work, while the 2007 observation (280 GHz) only employed seven antennae. As a result, despite 
its higher frequency, the 2007 data has a lower angular resolution than the 2008 data (as indi- 
cated by their synthesized beam sizes, in Table 3). We therefore used the 2008 data (frequency 
centered at 225 GHz, or 1.3 mm) to analyze the dust continuum emission. The continuum was 
averaged from the line-free channels and then subtracted from the side-band spectrum. The 
continuum data of the two sidebands were averaged on the (u,v) plane and then converted to 
the image domain. After Clean and Self- calibration, the 1.3 mm continuum image has an rms 
noise level (1 a) of 3.6 mjy beam" 1 (corresponding to a brightness temperature of It, = 0.0025 
K). 



http://irsa.ipac.caltech.edu/ 
7 http://www4.cadc-ccda.hia-iha.nrc-cnrc.gc.ca 

The Submillimeter Array is a joint project between the Smithsonian Astrophysical Observatory and the Academia Sinica Institute of Astron- 
omy and Astrophysics and is funded by the Smithsonian Institution and the Academia Sinica, see http: / /www.cfa.harvard.edu/sma/ 
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3 RESULTS 

3. 1 Dust continuum emission 

In Figure 4, we show the continuum emissions of G8.68 from infrared to (sub) millimeter wave- 
bands, including the IRAC 3.6, 4.5, and 8 |J.m emissions (RGB color image), the SMA 1.3 mm 
continuum emission (white contours), and the SCUBA 450 um continuum (blue dashed con- 
tours). 

Figure 4 is centered at the 1.3 mm continuum peak, the coordinates of which are RA.(J2000) = 
18 h 06 m 23.52 s , Decl.(2000)= -21°37'll". It is close in projection to the SMA phase tracking cen- 
ter (labeled with the red cross) . After deconvolution with the synthesized beam, the core has an 
angular size of 1 1" x 6" for the 4 a contour (0.24 x 0.13 pc at a distance of 4.5 kpc) . It is elongated 
in the north-south direction {PA = -10° for the major axis), reasonably coherent with the 4.5 
|j.m emission (green color). Since the 4.5 um emission traces the shock interaction between the 
outflow and the envelope gas, it is possible that the outflow and shocks are also affecting the 
dust distribution, causing its observed elongation. We did not find any evidence for multiple 
sub-cores either in our 1.3 mm continuum or any molecular lines (Figure 5). Therefore the gas- 
and-dust core should have a single compact morphology, and the fragmentation is not evident 
on our observational scale (0.05 to 0.5 pc). 

More diffused dust component can be revealed by the SCUBA 450 |j.m continuum emission. 
As shown in Figure 4, the 450 |J.m emission is more extended and less elongated than the 1.3 
mm emission. We use the average deconvolved FWHM (full width at half maximum) radius (r) 
to represent the extent of the continuum and molecular line emissions. Normally <r> can be 
measured from the 50 % contour level of the emission region. However, the 50 % contour (for 
the continuum and molecular lines) is often close to or even smaller than beam size. We thus 
suggest measuring the deconvolved radius from the 10 % contour, and adopt its 1/2 as the value 
of (r). For the continuum images, the 10 % contour is not specifically plotted in Figure 4, but 
close to the 14 a and 4 a contour level for the 1.3 mm and 450 |j.m emissions, respectively. We 
obtained (r) 450 ^, m = 0.23 pc and (r> L3mm = 0.08 pc. 

We also measured the integrated flux density F(X) of the dust core at wavelength A. In gen- 
eral, we use the 4 a emission level as the integration area for F(A). As an exception, for the IRAC 
data, we use the region of the 4.5 |xm emission (green color in Figure 4) to measure the inte- 
grated flux of all four bands, since the 4.5 |xm emission has a relatively clear boundary. The 5.8 
urn band shows a similar morphology with the 4.5 |xm, while the emissions at other two bands 
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are much fainter and cannot be well delineated. The IRAC stellar sources in the vicinity of the 
core are carefully excluded from the integration area. The derived F{\) are shown in Table 4. 



3.2 Molecular lines 

3.2.1 The Submillimeter Array 

Using the SMA, we have detected a number of molecular transitions of C 18 0, C 34 S, DCN, and 
CH 3 OH. Their beam-averaged spectra towards the 1.3 mm continuum peak and their velocity- 
integrated intensity images are shown in Figure 5. For the CH 3 OH, altogether we have detected 
1 1 rotational transitions. We selected five of them with largely different E u , and presented their 
images in Figure 5. The physical parameters of all the molecular transitions are shown in Table 
5. 

As shown in Figure 5, the spectra of the molecular tracers of high-density gas mostly show 
lines with single peak profiles. As exceptions, there are two CH3OH lines, II2 - IO3 {E u - 191 K) 
and 157,8 - 166,n {E u = 523 K) which show double peak profiles. However, since the remaining 
lines are all single-peaked, the two lines are more likely to be blended with other molecular 
transitions. The possible candidates for the blenders are NH 2 CN 14 2 - 13 2 (f=279. 35062 GHz, 
E u = 228 K) and HCCNC 29 - 28 (f=288.07346 GHz, E u = 207 K). We used two gaussian pro- 
files to fit the blended spectra, as plotted in dotted lines in Figure 5. For each spectrum, the 
peak velocity of the second component is well consistent with the anticipated velocity for the 
blenders. With the contamination excluded, these two CH3OH lines should also have single 
gaussian profiles. 

The C 34 S emission region shows an elongated morphology from the northeast to southwest 
{PA = 45°) as labeled in dashed line. The elongation agrees with the orientation of the CO out- 
flow and 4.5 urn shock emission (Figure 4). An elongated morphology towards northeast is also 
shown in the low excited CH 3 OH lines (i.e. E u = 33 K and 97 K). Therefore, both the CH 3 OH and 
C 34 S distributions should be affected by the outflow. The C 18 (2 - 1) also shows a non-regular 
morphology. However, it is biased to the south of the dust core, peaked at offset= (0",-2"), 
In addition, the C 18 shows a secondary clump in the southeast, peaked at offset= (10", -8"). 
This clump is not detected in C 34 S or CH3OH lines, indicating that it may be depleted in these 
species. The C 18 might trace cooler and less dense gas component, thus have a more extended 
feature than other dense-core species. For each molecular transition, we also measured the av- 
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erage deconvolved radius from the 10 % contour level (and adopted its 1/2 as the value of (r>). 
The results are shown in Table 5. 



3.2.2 The single dishes 

Figure 2 and 3 show the molecular lines detected from the single dishes. As shown in Figure 
3, prominent double-peak line profiles are observed in the 12 CO (1 - 0), (2 - 1), and also HCN 
(3-2). For both the 12 CO (1 - 0) and (2-1), the blueshifted emission peak is much stronger 
than the redshifted one, and the central absorption dip is well coincident with the C 18 line 
peak (Vj sr = 37 km s" 1 ). S uch blue asymmetric 12 CO lines sugge st the presence of infall motion 



towards the core center (Zhouet al 



1993 



Mardones et al 



19971) . For the physical explanation, 



when the infall occurs in the envelope which is cooler than the inner region, the gas in the front 
part would absorb the redshifted side of the line profile, whereas the gas in the rear part (behind 
the center) would increase the blueshifted emission because it is moving towards the observer. 
Besides the infall signature, the 12 CO lines also exhibit high-velocity emission wings extending 
to V[ sr =25 and 50 km s" 1 for the blue- and redshifted sides, respectively. This velocity range is 
comparable to the outflow velocities revealed by LI 1 (Figure 5 therein) . 

The two 13 CO lines also have a blueshifted emission peak (VJ sr = 35 km s _1 ) with respect 
to the C 18 0, suggesting that the 13 CO is also probably tracing the infall motion. However, be- 
cause the 13 CO lines are much less optically thick than the 12 CO, they exhibit no central dip, but 
instead show an emission shoulder that continuously declines towards the redshifted side. 

As shown in Figure 2a, we can see that the offset positions also exhibit self-absorbed pro- 
files (except the southeast one). However, compared to the central spectrum, their blue- and 
redshifted peaks have more similar intensities. As an extreme, the southeastern spectrum have 
a flattened top, with the double-peak feature almost disappeared. This indicates that the infall 
motion (along the line of sight) should have a decline towards those offset points. And their 
distance from the center (0.7 pc) can therefore be taken as a lower limit for the radius of the 
infalling region. 

As shown in Figure 3c, the HCN (3 - 2) has a double peak profile and high-velocity emis- 
sion wings extending to Vj sr = 29 and 53 km s _1 (above the noise level) for the blue and red 
wings, respectively. However, its double peaks have different asymmetry with the 12 CO lines, 
i.e., the red peak is slightly stronger than the blue one. The offset spectra of the HCN (Figure 
2b) have much lower signal-to-noise ratio (mainly due to their shorter integration time). How- 
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ever, they still evidently show double peak profiles, and have similar intensities for the blue and 
red peaks. Like in the case of the CO lines, the optically thin isotopic lines, i.e., DCN (3 - 2) and 
(4-3) are both single-peaked, hence the HCN (3-2) profile should originate from a self absorp- 
tion effect. Compared to the 12 CO profiles, the HCN spectra may r eflect different gas motions, 
including core expansion or/ and rotation ( Pavlyuchenkov et alj2008l) . In particular, as the most 
possible case, a cold and spherically expanding envelope would cause a prominent blueshifted 
self-absorption towards the center, while at the offset positions, the expansion should be on 
the plane of the sky, thus show a less blueshift due to the lower radial velocities. This scenario 
can reasonably explain the observed line profiles, but still needs a further examination with an 
improved angular resolution and spectral sensitivity. In Figure 2b, we also plotted the velocity- 
integrated map of HCN (3 - 2) (discrete gray scales) . We measured the average deconvolved 
radius <r> of HCN from its 50 % contour. As a result, it has <r> = 12" (and r/bf = !)• 

Another JCMT line, H 13 CO + (3-2) , has a regular gaussian profile, indicating that it may arise 
from the dense molecular core and is not evidently affected by the infall or outflow motion. 
We only have one-point observation for the H 13 CO + (3 - 2) at the continuum center, and in 
calculation for its column density (Section 4.4), we assume r/bf = 1- 



4 DISCUSSION 

4. 1 The physical properties of the dust core 

As shown in Figure 4, the 1.3 mm dust core does not coincide with any infrared stellar sources 
besides the extended 4.5 |xm shock emission. This indicates that the stellar emission from the 
core center is highly obscured by the dust. In the vicinity of the 1.3 mm dust core, a few stellar 
objects are shown in the IRAC RGB image (also labeled with the asterisks). All these objects are 
isolated from the 1.3 mm continuum emission, yet the three objects nearest to the continuum 
peak are likely to be embedded in the 450 u.m emission region. They might either be more 
evolved young stars in the same star forming region or irrelevant foreground stars. Despite this 
uncertainty, it is clear that these objects have no significant contribution to the dust continuum 
or molecular line emissions. We therefore make no further discussion for them. 

We can fit the Spectral Energy Distribution (SED) of the dust core from its flux densities at 
the Spitzer and JCMT/SCUBA wavebands. Assuming a gray body emissio n and a uniform dust 



temperature Id, the continuum flux density would be ISchnee et al. 



2007) 
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M core K v B v (r d ) 

F v = =r (1) 

gD 2 



where F v is the flux density at frequency v. M core is the total gas-and-dust mass of the core. 
g = 100 is commonly adopted gas-and-dust mass ratio. B v (T d ) is the Planck function at tem- 
perature Td- D = 4.5 kpc is the source distance. The dust opacity k v is assumed to have a 
pow er-law shape, i.e. k v = K 2 3q r T Hz(v/230GHz)^, with the reference value k 230 ghz = 0.9 cm 2 



g" 1 llOssenkopf & Henningl 19941) . The free parameters in the fit are M core , T&, and /}. We found 
that the emissions from 8 |j.m to 850 |J.m can be best fitted with two temperature components 
which have I'd = 20 K and 120 K respectively, and /3 = 2.1. The best-fit SED is shown in Figure 6. 

We did not include the IRAC 3.6, 4.5 or 5.8 urn emissions in our SED model. The 4.5 and 
5.8 |J.m emissions may largely come from the shocked emission thus are much stronger than 
the emissions at other two IRAC bands (Table 4). As for the 3.6 u.m emission, if being thermally 
excited, it may arise from some even hotter component which is much fainter and poorly con- 
strained by our current data. Therefore we also neglected the 3.6 u.m band. With the derived 
SED, the bolometric luminosity can be estimated using Lbol = 4ttD 2 f F v dv. As a result we have 
Lboi = 2.3 x 10 4 and 8 x 10 2 L e for the cold (20 K) and warm (120 K) components, respectively. 
Using Equation (1), we can also estimate the total mass of the two temperature components, 
which turn out to be 1.3 x 10 3 and 1.0 x 10~ 3 M Q for the 20 K and 120 K components, respec- 
tively. One can see that both the core mass and luminosity are dominated by the gas-and-dust 
component which is characterized by T d = 20 K. 

In Equation (1), by replacing the integrated flux density F v with the flux density at the con- 
tinuum peak (0.32 Jy beam" 1 ), and then dividing the obtained mass with the beam area and 
the average molecular mass (1.4 times the molecular mass of H 2 ), we can derive the H 2 column 
density N (H 2 ) towards the continuum peak. And then, assuming that the core is approximately 
spherical, we can derive the volume number density using n(H 2 ) = A/"(H 2 )/2(r). The physical 
parameters of the dust core are presented in Table 4. 

By extrapolating the best-fit SED curve, we can get a flux density of 2.9 Jy at A = 1.3 mm. 
Compared to this value, the SMA observation has recovered 35% of the 1.3 mm continuum 
emission. Adopting T d = 20 K, we also estimated the physical parameters from the SMA 1.3 mm 
continuum, which are presented in Table 4. Based on the continuum observations, we suggests 
that the gas-and-dust core in G8.68 should consist of a dense inner region (characterized by the 
1.3 mm emission), and a more extended envelope (traced by the 450 urn emission). 
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4.2 The CH3OH rotational temperature 

The molecular gas temperature can also be estimated from the CH 3 OH lines using the rotation 
diagram. The methanol lines all have linewidths of several km s~ : , with none of them showing 
abnormally high intensities, therefore the CH3OH lines are unlikely to have maser excitations. 
Assuming optically thin, th e column de nsity of the upper-level N u can be derived from the 



2005) 



integrated line intensity using (Tielens 
87ikv 2 , r 

N ° = hc^} T » dvl ™ (2) 

where It, is the observed brightness temperature. A u \ is the Einstein coefficient in s" 1 . r/bf is 

the beam filling factor. All other constants take their usual values in SI units. Although for an 

emission line, rj^f may vary at different velocities, we approximate it to be a single value as the 

ratio between the integrated emission region and the beam area, i.e., rj^f ^ 7i{r) 2 l 'Abeam (used 

when the emission region is smaller than the beam size, otherwise r/^f = 1). r^bf is estimated for 

each transition and the derived values are presented in Table 5 (Column 11). 

Assuming a Local Thermal Equilibrium (LTE, i.e. energy levels are populated according to a 

Boltzmann distribution characterized by a single temperature), the relation between the total 

column density Ny and N u is 

N u N T E u 

— = 7^r^exp(-— ^) (3) 

gu Qdrot) kT Iot 

and its logarithmic form is 

N u N T E u 

ln( — ) = ln(— — -) - — — (4) 

gu Qdrot) kT mt 

where g u and E u are the degeneracy and the excitation energy of the upper level, respectively 
and Q(r r nt) is the partition fun ction. For CH3OH, a good approximation is Q(T Iot ) - 1.2327 x 
T}^ jTownes & Schawlowlll955l) . 



rot 



The rotation diagram for the CH 3 OH lines is shown in Figure 7. A linear least-square fit to 
the data points results in T mt = 130 + 10 K and N T = (5.3 ± 0.6) x 10 15 cm" 2 . 

In the calculation, in order to correct for the optical depth effect, one should multiply N u /g u 
with a correction factor C V = r/(l-e~ T ) andfit the rotational temperature iteratively. The optical 



depth is estimated using JRemijan et alj|2004l Equation (3) therein, slightly reformed) 



c 3 V4ln2 hv 

r -^w^ NMexpi ^ ) - li (5) 

Among all the CH3OH lines, the (8-1 - 7 ) transition has the highest optical depth (t = 0.059). 
The other lines are even more optically thin. To take into account the temperature uncertainty, 
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we also estimated the optical depth assuming T mt = 20 K which is a lower limit as suggested by 
the SED fitting. At 20 K, the optical depths become ~ 8 times larger than the values at T rot = 130 
K. The derived optical depths are listed in Table 5, and the column densities and abundances 
are listed in Table 6. 

The rotational temperature of T I0t = 130 K is close to the SED temperature of the warm dust 
component (Id = 120 K). Therefore it is possible that the CH3OH emissions are mainly from the 
region associated with the warm dust. Moreover, since the cold component (Id = 20 K) is more 
massive than the warm one for orders of magnitude, the CH3OH may have a severe depletion 
in the region for the cold dust component. However, it is also possible that the dust and gas 
are thermally decoupled, thus exhibit different temperatures. The collisional excitations of the 
molecular gas can be particularly enhanced by the shocks (especially along the outflow direc- 
tion), thereby showing a high value of r rot . The dust temperature T&, in comparison, may still 
be largely dominated by the stellar heating thus has a much lower value. 

4.3 The CO emission and the infall motion 

As shown in Section 3.2.2, both the infall and outflow signatures are detected in the 12 CO (2 - 1) 
and (1 - 0) lines. In this paper we mainly discuss the infall properties based on the 12 CO (2 - 1). 
We first make attempt to separate the different components from the observed spectrum, then 
estimate the infall rate. 

Following the procedure of ! 



Purcell et al. 



2006), we used a broad gaussian profile to fit the 



outflow wings (velocity range of V < 32 and V > 42 km s _1 ), and then subtracted it from the 
spectrum. The residual line profile (green line in Figure 8) should mainly represent the emis- 
sion from the dense molecular core. One can then mask the velocity range possibly affected by 
the infall motion (34 to 43 km s _1 ), and make a gaussian fit to the spectrum outside this veloc- 
ity range. The fitted spectrum is speculated to roughly represent the molecular core emission 
unaffected by the infall signature. However, for the 12 CO lines, due to its large optical depth, we 
cannot directly apply a Gaussian fit to the spectrum. Instead, one should model the spectrum 
using the radiation transfer function. In this case, the line profile can be expressed as 

T mh {V) = [T mh , -KT CMB )][l-e- T[V) ] (6) 

J{T) = r /[exp(r /r)-l] is the Planck- corrected brightness temperature, and T = hvlk. Xcmb = 
2.7 Kis the temperature of the cosmic background. At the frequency of CO (2-1) (230 GHz), we 
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have /(Icmb) = 0.2 K. Compared to the intensity of the CO emission, the contribution from the 
cosmic background can be almost neglected. 

We also assume the dense molecular core to have a uniform gas distribution along the line 
of sight, with central velocity V , velocity dispersion a and peak optical depth t . Then the 
optical depth is 

T(tO = T exp[- 2 J J ] (7) 



where a is related to the (intrinsic) line width A V by a = AV7\/8ln2. In an optically thick case, 
the line emission could be largely saturated. To is thus poorly constrained by the observed 
spectrum. Howeve r, it can be estimated from comparison to the CO isotopologues following 



Garden et al 



(1991 



, Equation (4) therein) . Since the CO is also affected by the infal l motion, we 



used the C 18 (2 - 1) instead. Assuming an abundance ratio of [ 12 CO/C 18 0] = 490 llGardenet al. 



19911) . the equation will be 

r mb ,o( 12 CO) l-exp[-T ( 12 CO)] l-exp[-T ( 12 CO)] 



(8) 



r mb ,o(C 18 0) l-exp[-T (C 18 O)] l-exp[-T ( 12 CO)/490] 

To fit the line profile, we first take an arbitrary, but reasonable value of To, and then fit the 
line profile by adjusting the values of r m b,o, V and AV in Equation (6) and (7). The best-fit 
T m b,o is then used to estimate t again using Equation (8). The final best fit can be reached 
after two or three iterations. Eventually, we have r m b,o = 32 K, AV = 5.5 km s~\ Vb = 37 km 
s"\ and To( 12 CO) = 68. In Figure 8, the best fit spectrum is shown in dashed line. And a sum of 
dense-core and outflow components is shown in red dot-dashed line. The output spectrum has 
an apparent line width of 8.5 km s _1 which is indeed much broader than the intrinsic A V. We fit 
the 12 CO (1 - 0) using the same method. All their line parameters are listed in Table 5 (Column 
5 to 8). 



The infall rate is estimated using IIKlaassen & Wilson 



2007) 



M in f=-7rn(H 2 )iJ.miir gm Vi n (9) 

wherein r gm is geometric mean radius of the core, n(H 2 ) is the ambient source density, and 
Vi n is the typical infall velocity. I n calculation we es timated V\ n from the outflow-subtracted 



line profile using Equation (9) in 



Myers et al. 



(1996). As a result we have Vj n = 0.45 km s" 



In addition, we assume that the more diffused gas traced by the 450 urn emission which has 
rc (H2) = 0.8 x 10 6 cm -3 , is collapsing towards the dense inner region characterized by the 1.3 
mm continuum (Figure 4), thus we have r gm = (r) 1.3mm = 0.08 pc. With these assumptions, we 
derived an infall rate of M in f = 7.0 x 10" 4 M yr" 1 . 
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As seen in Equation (9), the derived infall rate is sensitive to the adoption of r gm , and our 
currently adopted r gm is relatively conservative. Adopting r gm = (r) 45 o(j.m = 0.23 pc, we would 
have M in f = 5 x 10" 3 M Q yr" 1 . However, we note that such a large-scale estimate may deviate 
from the small-scale infall rate. To resemble the mass infall onto the central stars, it may be 
more reasonable to adopt the first value (r gm = 0.08 pc). With the obtained infall rate, we then 
estimate the accretion luminosity, using L acc = GM,M in f/fi*, and assuming a mass-radius re- 
lation of R*IR Q = (M* /M ) 08 . As a result, we have L acc = (4 + 2) x 10 4 L e . The uncertainty in 
L acc corresponds to a stellar mass varying between 10 an 100 M . It is likely that the bolometric 
luminosity of the dust core (2.3 x 10 4 L Q , see Section 4.1) should have a major energy supply 
from the accretion process. 

Considering the existence of the outflow, there should be a strong interaction between the 
infall and the ou tflow. And the interaction may be responsible for the 4.5 |xm shock emission. 



Chen et al 



1 20101) have performed an HCO + (1-0) survey towards the Extended Green Objects 
(EGOs), i.e., the massive YSO candidates with the 4.5 um shock emissions. They found that 
nearly one third of the sample (29 out of 69 sources) exhibit a significant blue asymmetry. While 
in an HCO + (1-0) survey towards 8 2 massiye YS Os selected from the methanol masers, only 12 



sources have infall signatures flPurceil et al.ll2006l) . Comparing these results, it is likely that the 



shocks are prone to take place in the YSOs with infall motions. This is theoretically expected, 
since compared to an interaction between the outflow and quiescent gas, an outflow-infall in- 
teraction would more efficiently convert the kinetic energy into heat and radiation. 



4.4 The molecular abundances and DCN overabundance 

The total column density AT T (X) and abundance /(X) = N T (X)/N(H 2 ) of the C 18 0, HCN, DCN, 
H 13 CO + , and C 34 S are calculated from their emission lines at the continuum emission peak 
using equation (2) and (3). And a correction for the optical depth is done using Equation (4). 
To derive the HCN and H 13 CO + abundances, we used the AT(H 2 ) value for the SCUBA 450 u.m 
continuum (Table 3). While for the SMA lines, we adopted iV(H2) from the 1.3 mm continuum 
(0.95 x 10 24 cm" 2 ). We also note that A/t(H 13 CO + ) may be underestimated due to the beam 
dilution thus should be regarded as a lower limit. To take into account the temperature uncer- 
tainty, we also estimated AT T (X) at the lower limit of r rot = 20 K (suggested by the SED fitting). 
The Nj(X) and /(X) values are shown in Table 6. 

In calculation of the HCN, its line profile should be corrected for the self absorption. We 
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modeled its original line profile using the same method for the 12 CO (2 - 1) (Section 4.2). How- 
ever, since the abundance ratio between DCN and HCN is much more uncertain than [C 18 0/ 12 CO], 
we cannot use DCN to reliably determine the optical depth of HCN (3 - 2) . Nevertheless, we ex- 
pect the HCN (3 - 2) to have a low optical depth due to two reasons. Firstly, since the HCN (3 - 2) 
likely traces denser and hotter gas than the 12 CO (2 - 1 ), if the HCN (3 - 2) has a very large optical 
depth, it should have a comparable intensity with the 12 CO (2-1). Nevertheless, even after the 
self- absorption correction, the HCN (3 - 2) is still much weaker than the 12 CO (2-1). Second, 
with an apparent line width (AV = 6.2 km s~ 2 ) being close to AV of the C 18 and CH3OH lines 
(as shown in Table 5), the optical-depth broadening should be insignificant for the HCN (3-2). 
We therefore directly calculate iVr(HCN), and then estimate the optical depth using Equation 
(5). As a result, we found r = 0.78 at T Iot = 20 K and 0.10 at 130 K. This result is consistent with 
our expectation. However, to more accurately determine the HCN optical depth, one should 



consider to observe some other isotopologues such as HC N iHatchell et al 



1998) 



From the derived abundances, we can get a relative abundance ratio between DCN and 
HCN which is [DCN/HCN]^ 0.07. The values derived at the two t emperature limits are similar 



to each other (Table 6). Compared to the cosmic [D/H] ratio (10~ 5 . lLinslcvll998l] .the [DCN/HCN] 



in G8.68 implies a deuteration for orders of magnitudes. The [DCN/ HCN] in G8.68 is also much 



Hatchell et al 



1998). How- 



higher that the values detected in hot molecular cores (10 4 to 10" 

ever, it is much more comparable to t he abundance rat io of [N 2 D + /N 2 H + ] detected in high 



mass YSOs in the infrared dark clouds (Chenetal. 



detected in a number low-mass YSOs (Roberts et al 



2011). Overabundant DCN was previously 



2002), while in high-mass star-forming re- 



gions, the DCN is not frequently detected. 

The overabundant DCN in G8.68 may originate from a high level of deuterium fractionation 
in the previous cold pre-stellar phase. In highly deuterated gas (abundant in H2D + , CH2D + , 
C2HD + etc.), DCN can be produced via D-H s ubstitution of the HCN, or through more com- 



plex reaction pathways flAlbertsson et al.ll201ll. Reacti on (17) to (21) therein). Finally the DCN 



molecules would mostly reside on the grain mantles IIHatchell et al 



1998; 



Roberts et al. 



2002) 



During the star formation, the DCN can be released into the gas-phase again. However, once 
the temperature slightly increases, the gas-phase DCN can be easily destroyed via reactions 
such as 



H + DCN^HCN + D Oiarnley et al 



1992; 



Roberts et al. 



2002, Figure 5 therein). In this sense, 



to maintain the DCN overabundance in the gas, two conditions may have to be satisfied. First, 
there should be a high-level deuterium fractionation in the previous dark-cloud phase. Second, 
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in order to compensate the chemical destruction due to the stellar heating, a rapid sublimation 
for the dust grains should be necessary. Again, the outflow and shocks may have a potential 
contribution to this process. However, unlike the C 34 S and CH 3 OH, the DCN emission has a 
compact and spherical morphology which is not evidently coherent with the outflow. Therefore 
it is possible that the DCN enhancement is less affected by the shocks and/ or more sensitive to 
the stellar heating. A higher sensitivity and spatial resolution may help better reveal the DCN 
morphology and determine whether it is associated with the outflow. 

As another po ssibility, the DCN can also be synthesized in the recent gas phase chemistry. 



Parise et al 



2009) show that the gas-phase reactions may sufficiently account for the enhanced 
D-H ratio in the molecular gas in the Orion Bar PDR, which has a [DCN/HCN] ratio of 10 -2 , 
comparable with that in G8.68. However, the gas-phase enhancement may have to proceed 
in an environment with stable lukewarm heating. This condition may hardly be satisfied in 
regions with rapid evolution of the high-mass stars. Therefore the grain mantle sublimation 
may still be the major process for the DCN enhancement in G8.68. In the future study, one 
can compare other chemical products from the grain sublimation and gas-phase chemistry in 
order to determine the relative importance of these two processes. 

The C 34 S abundance in G8.68 i s much lower than the average C 34 S abundance in the UC 



Olmi & Cesaroni 



Frereking et al 



1999). /(CO) is also much smaller than the typical ISM value 



19821) . Compared to the ISM abundance, it has a depletion factor of 



Hn regions ( 
(1.7 xlO -7 , 
/ D (C 18 0) = 5±2. 



5 SUMMARY 

We have investigated the dust continuum and molecular line emissions towards the high mass 
YSO G8. 68-0.37. We have revealed a dense compact gas-and-dust core in the SMA 1.3 mm con- 
tinuum emission, and its more extended envelope in the SCUBA 450 urn emission. At our an- 
gular resolution (spatial scale > 0.05 pc), there is no evident fragmentation structures. We find 
that an SED with at least two temperature components is necessary to account for the dust con- 
tinuum emissions from mid-IR to sub millimeter wavelengths. The best-fit temperatures for the 
two components are Id = 20 K and 120 K. The core mass and luminosity are mainly contributed 
by the cold component (Id = 20 K). 

Prominent infall signatures and outflow wings are detected in both 12 CO (1 - 0) and (2 - 1) 
lines. We separated the outflow and dense-core components and measured their line param- 
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eters. The 12 C0 (2 - 1) yields an infall velocity of 0.45 km s _1 . Assuming that the extended en- 
velope is collapsing towards the inner dense region, we can derive an infall rate of 7 x 10~ 4 M Q 
year" 1 . It is possible that the 4.5 luti shock emission is largely enhanced by a strong interac- 
tion between the infall and outflow motions. In addition, we also suggest that the infall motion 
may be important for suppressing the stellar emissions thereby protecting the DCN and other 
fragile species. 

We estimated a rotational temperature of 130 K from the CH3OH lines. We derived the abun- 
dances of the molecular species from their spectra, and in particular, we found a high abun- 
dance ratio of [DCN/HCN]= 0.07. The over abundant DCN may originate from a high-level of 
deuterium fractionation in the previous pre-protostellar phase, as well as the recent grain man- 
tle sublimation and/or gas-phase chemistry. More details in this chemical process are still to be 
further investigated. 
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Table 1. General information of the single dish observations. 



Instrument 


PMO 






CSO 






JCMT 


Obs. date 


Dec 2009 






May 2011 






Aug 2005 


Beam size 


56" 






30" 






21" 


Vmb 


0.62 






0.698 






0.63 


Pointing center 


R.A.= 18:06:22.87 


R.A. 


= 18:06: 


23.5 


R.A. 


= 18:06:23.46 




Dec.= -21:37: 


:20.7 


Dec. 


= -21:37 


:10.7 


Dec.; 


= -21:37:09.64 



Table 2. observational parameters for the molecular lines from the single dishes. 



Transition 


Instrument 


Atmosphere 


Band width 


AVres 


rms noise 






opacity 


(MHz) 


(kms -1 ) 


per channel (K) 


12 CO (2 - 1) 


CSO 


0.167 


500 


0.079 


0.2 


13 CO (2 - 1) 


CSO 


0.152 


500 


0.083 


0.2 


C 18 0(2-l) 


CSO 


0.149 


500 


0.083 


0.2 


12 CO(1-0) 


PMO 


0.015 


145 


0.370 


0.1 


13 CO(1-0) 


PMO 


0.015 


43 


0.115 


0.1 


C 18 O(l-0) 


PMO 


0.015 


43 


0.115 


0.2 


HCN (3 - 2) 


ICMT 


0.111 


160 


0.088 


0.5 


H 13 CO + (3-2) 


ICMT 


0.067 


160 


0.090 


0.3 



Table 3. Observational parameters of the SMA. 



Epoche 



Frequency bands (GHz) 
LSB, USB 



Bandpass Flux Phase & synthesized beam size rms noise 

Calibrator Calibrator Calibrator (arcsec) per channel (K) ° 



2007 (279.4,281.4), (289.4,291.4) 

2008 (219.5,221.5), (229.5,231.5) 

2009 (217.5,219.5), (227.5,229.5) 



3c273 


Neptune 


1733-130,1911-201 


7.0 x 5.8 


0.110 


3c454.3 


Neptune 


1733,1911 


4.8 x 3.6 


0.017 


3c273 


Uranus 


1733,1911 


6.8 x 3.6 


0.150 



a. For the unit conversion, 1 K=0.367, 1.43 and 1.04 ly beam for the data in 2007, 2008, and 2009 respectively. 
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Table 4. Physical parameters of the dust core. 



Parameter 


Value 


Unit 


F(3.6um) 


13 + 0,1 


mjy 


F(4.5um) 


78±2 




F(5.8um) 


68±2 




F(8.0um) 


10 + 1 




F (24 urn) 


657 + 40 




F (70 urn) 


196 + 15 


Jy 


F (450 urn) 


144 + 15 




F(850|xm) 


15 + 4 




F(1.3mm) fl 


0.65 + 0.03 




F(1.3mm) b 


2.9 






(450 u-m / 1.3 mm) c 




<r> 


0.23 + 0.05 / 0.08 + 0.02 


pc 


j^core 


1.5 + 0.2/0.30 + 0.01 


1O 3 M 


iV(H 2 ) 


1.2 + 0.2/ 0.95 ±0.03 


10 24 cm -2 


«(H 2 ) 


0.8 + 0.1 / 3.8 ±0.3 


10 6 cm" 3 



Note. To measure the integrated flux density of the dust core, we 
use the 4 a emission level as the integration area (aperture for the 
photometry). As an exception, we use the 4 a level of the 4.5 u.m 
emission as the area for all four IRAC bands, which roughly equals to 
the emission region with green color in Figure 4. The nearby point 
sources carefully excluded from this aperture. 

a. The flux density of the SMA continuum observation. 

b. The flux density extrapolated from the SED fitting (Figure 6). 

c. For the last 4 parameters, the first and second values are derived 
from the 450 um and 1.3 mm continuum data, respectively. 
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Table 5. Observed parameters of the molecular lines. 



Molecule 


Transition 


Rest frequency 


E u 


Vlsr 


^b.peak 


A^FWHM 


fT h dv 


T a 


(r) b 


<i 








(GHz) 


(K) 


(kms" 1 ) 


(Kelvin) 


(kms" 1 ) 


(Kkms" 1 ) 




(arcsec) 




(1) 


(2) 




(3) 


(4) 


(5) 


(6) 


(7) 


(9) 


(8) 


(10) 


(11) 


12 CO(core) 


1-0 




115.27120 


5.5 


37.0 


26.0(1.5) 


5.0(0.5) 


134(15) 


75 


- 


- 


12 CO(outflow) 


1-0 








37.3 


4.0(1.5) 


9.2(1.0) 


37(5) 


0.15 


- 


- 


12 CO(core) 


2-1 




230.53800 


17 


37.0 


32.0(0.5) 


5.5(0.5) 


240(30) 


68 


- 


- 


12 CO(outflow) 


2-1 








37.0 


5.0(0.5) 


9.0(0.5) 


45(6) 


0.15 


- 


- 


HCN(core) 


3-2 




265.88643 


26 


38.5 


12.0(0.2) 


6.2(0.5) 


73(9) 


0.78/0.10 


12 


1 


HCN(outflow) 


3-2 








39.5 


0.8(0.2) 


15.0(1.0) 


12(5) 


0.053/0.007 


- 


- 


13 CO 


1-0 




110.20135 


5.3 


35.0 


12.4(0.5) 


6.8(0.1) 


85(10) 


0.85 


- 


- 


13 CO 


2-1 




220.39968 


16 


35.0 


8.5(0.5) 


7.5(0.5) 


60(15) 


0.70 


- 


- 


C 18 


1-0 




109.78217 


5.3 


37.0 


4.5(0.5) 


5.6(0.4) 


23(4) 


0.17 


- 


- 


C ls O(CSO) 


2-1 




219.56036 


16 


38.0 


6.0(0.3) 


5.5(0.4) 


34(3) 


0.140 


- 


- 


C 18 0(SMA) 


2-1 




219.56036 


16 


38.0 


5.2(0.02) 


4.8(0.3) 


30(3) 


0.423/0.065 


3.4 


1 


CH 3 OH 


80 -7i 


220.07849 


97 


37.5 


1.6(0.22) 


5.0(1.3) 


8(2) 


0.209/0.028 


3.1 


1 


CH3OH 


8-1- 


7 


229.75880 


89 


37.0 


3.42(0.15) 


6.0(0.3) 


18.0(1.5) 


0.453/0.059 


2.8 


1 


CH3OH 


3-2- 


•4-1 


230.02706 


39 


38.5 


0.9(0.03) 


6.0(0.4) 


5.8(0.5) 


0.036/0.004 


3.2 


1 


CH3OH 


10 2 - 


9 3 


231.28110 


165 


39.0 


0.53(0.06) 


6.0(0.9) 


3.4(0.5) 


0.087/0.011 


2.2 


0.95 


CH3OH 


9-1- 


80 


278.30451 


110 


38.0 


1.39(0.06) 


5.4(0.4) 


7.2(0.8) 


0.245/0.031 


2.9 


0.89 


CH3OH 


2-2" 


•3-1 


278.34226 


33 


39.0 


0.24(0.02) 


5.5(0.4) 


1.5(0.2) 


0.045/0.006 


3.2 


1 


CH3OH 


21-2 


-20_3 


278.48023 


563 


40.0 


0.09(0.02) 


4.0(0.5) 


0.4(0.05) 


0.024/0.002 


2.2 


0.48 


CH3OH 


157,8 


- 166,11 


288.07677 


523 


37.5 


0.24(0.02) 


5.5(0.3) 


1.2(0.2) 


0.102/0.012 


2.3 


0.52 


CH3OH 


144- 


153 


278.59906 


340 


39.5 


0.29(0.02) 


6.3(1.3) 


1.8(0.2) 


0.073/0.009 


2.5 


0.62 


CH3OH 


11 2 - 


10 3 


279.35191 


191 


39.5 


0.28(0.02) 


6.5(1.9) 


2.3(0.3) 


0.081/0.009 


2.8 


0.78 


CH3OH 


43 -52 


288.70557 


71 


39.5 


0.28(0.02) 


6.5(1.5) 


1.8(0.3) 


0.068/0.008 


2.8 


0.77 


DCN 


3-2 




217.23854 


81 


39.0 


2.1(0.4) 


3.5(0.2) 


7.5(1) 


0.132/0.018 


3.1 


0.92 


DCN 


4-3 




289.64492 


35 


39.0 


1.0(0.04) 


3.5(0.2) 


5.0(0.4) 


0.072/0.009 


3.3 


1 


C 34 S 


6-5 




289.20907 


38 


39.0 


1.4(0.02) 


4.5(0.5) 


6.5(0.7) 


0.123/0.016 


3.3 


1 


H 13 CO + 


3-2 




260.25534 


25 


38.0 


4.7(0.3) 


4.4(0.3) 


20(2) 


0.260/0.033 


- 


- 



Note. The 12 CO, 13 CO lines, and C 18 lines are from the PMO and CSO observations (see Table 2). The C 18 (2 - 1) line from the SMA 
observation is also presented. The HCN and H CO* lines are observed with the JCMT. For the double-peaked lines, including HCN (3—2), 
CO (2 - 1) and (1-0), and two CH3OH lines (£ u = 191 K and E u = 523 K), the parameters are measured from the fitted spectra, while for 
the single-peak lines, we directly measured the observed spectra. 

a. The optical depth at the line peak. For the transitions with two values, the first and second one are the results for r ro t = 20 K and 114 
K, respectively (see Section 4.4). While for the CO (1 - 0) and (2 - 1), the optical depths are calculated from comparing their isotopic lines 
(Section 4.3). 

b. The effective radius of the emission region, measured from the deconvolved average radius of the 10% contour region (1/2 times the 
value). An exception is the HCN (3-2), for which we directly measured the 50 % contour. The average uncertainty level is ~ 2 arcsec. For 
G8.68 at D = 4.5 kpc, 1 arcsec= 0.02 pc. 

c. The beam filling factor, calculated from the ratio between the area of the deconvolved emission region iji(r) ) and the beam size. 

Table 6. Collum density and abundance of the molecular species. 





Trot 1 


= 20K a 




Trot : 


= 130K 


Molecule 

(X) 












iV T (X) (cm~ 2 ) 


fOQ 




N T (X) (cm -2 ) 


fOQ 


c 18 o 


(2.2±0.3)xl0 16 


(2.3 ± 0.3) x 10" 


-8 


(5.4 + 0.6) x 10 16 


(5.6 + 0.7) xlO -8 


CH 3 0H fe 


(1.8±0.2)xl0 15 


(2.0 ± 0.2) x 10" 


-9 


(5.3 + 0.6) x 10 15 


(5.8 + 0.6) xl0~ 9 


C 34 S 


(2.3±0.3)xl0 13 


(2.5 ± 0.4) x 10" 


-11 


(2.6 + 0.3) x 10 13 


(2.9 + 0.4) xl0~ n 


H 13 CO + 


(9.5±1.0)xl0 12 


(7.9 ± 0.7) x 10" 


-12 


(1.7 + 0.2) x 10 13 


(1.4 + 0.1) xlO -11 


HCN 


(5.6 + 0.5) xlO 14 


(4.6 ± 0.3) x 10" 


-10 


(1.1 ±0.1) x 10 15 


(8.9 + 0.8) xlO -10 


DCN 


(3.1±0.1)xl0 13 


(3.4 + 0.2) x 10" 


-11 


(5.6 + 0.3) x 10 13 


(6.2 + 0.5) xlO -11 


[DCN/HCN] C 


- 


0.07 + 0.01 




- 


0.07 + 0.01 



a. A lower limit as suggested by the dust continuum SED. 

b. At T I0t = 20 K, N T (CH3OH) is derived from 3_ 2 - 4_i line. 

c. Abundance ratio between DCN and HCN. 
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Figure 1. The beam size and pointing center of each instrument. The cross symbol marks the center of the 1.3 mm dust core (Figure 4) which 
is coincident with the CSO pointing center. The white ellipse is the synthesized beam of the SMA in the 2008 observation, and the gray filled 
circle is the primary beam. The fCMT beam size is for the frequency of 289 GHz. 
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Figure 2. (a) Grid spectra of CO (2 - 1) observed from the CSO. The red cross labels the pointing center of each spectrum. The gray contours 
are the SCUBA 450 yim emission (specified in Figure 4) . The green dashed line represents the beam size, (b) Grid spectra of HCN (3 - 2) observed 
from the JCMT The red cross labels the pointing center of the each spectrum. The intensity map (gray scales) is made from interpolating the 
line intensity at each point. The integration for the spectra is from 25 to 50 km s . The gray-scale levels are from 30 % to 90 % of the maximum 
intensity (46.8 K km s _1 ). The thick contour is the 50 % level. The blue dashed circle is the beam size. 
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Figure 3. Molecular lines observed from the PMO, CSO, and 1CMT, which are shown in left, middle, and right panels, respectively. The observing 
centers and the beam size of each telescope are shown in Figure 1. 
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Figure 4. Continuum emissions detected towards G8.68 from infrared to millimeter wavelengths. The image is centered at the emission peak of 
the 1.3 mm continuum, the coordinates of which are RA.(J2000)=18 h 06 m 23.5 s and Decl.(J2000)=-21°37'l0.7". The white contours are the 1.3 
mm emission observed from the SMA. The contour levels are -4, 4, 14, 24... 104 a (0.003 Jy bea m ), The -4 a co ntour is due to the insufficient 
(u,v) coverage and is plotted in dotted line. The square denotes the strongest CH3OH maser iWalsh et al. 1998). The dashed contours are the 
JCMT/SCUBA 450 urn emission. The levels are 4, 8, 12... 36 o (1.2 Jy beam -1 ). The IRAC 3.6 (blue), 4.5 (green) and 8.0 (red) urn images are 
shown together in RGB colors (also seen in Figure 1 and 2 in LI 1) . The red cross labels the SMA phase center. The synthesized beam of the SMA 
(white ellipse) and SCUBA (blue circle) beam are plotted in the right corner. 
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Figure 5. Molecular lines and integrated images observed from the SMA. For each transition, the contours are 10, 20... 90 percent of the peak 
intensity. The integration range is (37, 43) km s for all the transitions except the two blended CH3OH lines. For these two lines the integration 
range is (37,40) km s as to eliminate the contamination. The dashed line in the C S labels the orientation of the outflow in Lll. The DCN 
(3 - 2) line is shifted above the zero level for 2 K. The negative contours due to the missing flux are omitted to more clearly show the emission 
features. The gray-scale image in each panel is the SMA 1.3 mm continuum. 
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Figure 6. The spectra energy distribution of the dust core. The black squares with error bars are the data points (the SMA 1.3 mm flux density 
is marked with the diamond) . The black line is the fitted SED curve. The red dashed lines are the SEDs of the two temperature components, 
with r d = 20 K and 120 K, respectively. 
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Figure 7. The rotation diagram of the CH3OH lines. The black squares with error bars are the data points. The red line is the least-square fit. 
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Figure 8. The two-component fit to the CO (2 - 1) spectrum towards the center. The meaning of each line type is labeled in the legend. 
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